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ABSTRACT 

B[e] supergiants (B[e]SGs) are evolved massive stars in a short-lived transition phase. During this 
phase, these objects eject large amounts of material, which accumulates in a circumstellar disk-like 
structure. The expelled material is typically dense and cool, providing the cradle for molecule and dust 
condensation and for a rich, ongoing chemistry. Very little is known about the chemical composition of 
these disks, beyond the emission from dust and CO revolving around the star on Keplerian orbits. As 
massive stars preserve an oxygen-rich surface composition throughout their life, other oxygen-based 
molecules can be expected to form. As SiO is the second most stable oxygen compound, we initiated 
an observing campaign to search for first-overtone SiO emission bands. We obtained high-resolution 
near-infrared L-band spectra for a sample of Galactic B[e]SGs with reported GO band emission. We 
clearly detect emission from the SiO first-overtone bands in CPD-52 9243 and indications for faint 
emission in HD 62 623, HD 327 083, and GPD-57 2874. From model fits, we find that in all these stars 
the SiO bands are rotationally broadened with a velocity lower than observed in the GO band forming 
regions, suggesting that SiO forms at larger distances from the star. Hence, searching for and analyzing 
these bands is crucial for studying the structure and kinematics of circumstellar disks, because they 
trace complementary regions to the GO band formation zone. Moreover, since SiO molecules are the 
building blocks for silicate dust, their study might provide insight in the early stage of dust formation. 
Subject headings: circumstellar matter — infrared: stars — stars: early-type — stars: massive — 
supergiants 


1. INTRODUGTION 

Massive stars, when evolving off the main sequence, 
undergo phases of strong mass-loss, often resulting in 
the formation of shells, disks or rings of circumstellar 
material. One class of evolved massive stars, B[e] super¬ 
giants (B[e]SGs), are early-type emission line stars with 
high-density gaseous and dusty circumst ellar disks of yet 
unknown origin (see de Wit et al. 2014, for a recent re¬ 
view) . The gaseous inner disk is traced by the emission 
of low-ionized metal lin es from both permitt ed and for¬ 
bidden transitions (e.g., Zickgraf et al. 1985). Of these, 
the lines of [Gall] and [OIJ are particularly useful. While 
their intensities probe regions of high and medium densi¬ 
ties, where the transition between the two regimes might 
be settled at an electron density of ~ 10^cm“^, their 
line profiles contain the information on their kinemat¬ 


ics (Kraus et al. 2007 2010 
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emission feature s in optical spectra of several B[e]SGs ha s 
been suggested ( Zickgraf et al.|[T9^ Torres et al.||2012 ) 


but still lacks conhrmation. Besides atomic and molec- 
ular gas, these stars are also surrounded by warm cir¬ 
cumstellar dust, as is obvious from their strong infrared 
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group of early-type blue supergiants (BSGs) with ring 
nebulae such as, e.g., Sher25 and SBW1. Optical im¬ 
ages of these t wo stars disp^y equ atorial rings and bipo¬ 
lar lobes, and Smith et al. (20071 propose that as soon 
as the rings or disk-like structures seen around B[e]SGs 
have expanded and cooled they would look exactly like 
these objects. Support for such a possible evolution¬ 
ary link comes from the fact that, like B[e]SGs, neither 
Sher 25 nor SBW 1 has evolved yet through a red su¬ 
pergiant phase, indicating that in both cases, the ejec¬ 
tion of the material forming the disks or rings must 
have happened during their BSG phase. Such a sce¬ 
nario has a much wider impact: as the BSGs with rings 
closely resemble the progenitor star of the supernova SN 
1987A, B[e]SGs might be regarded as the most plausible 
candidates for supernova explosions of the SN 1987A- 
type. Notably, two B[e]SGs were recently suggested to be 
viable SN 1987A-type progenitor can didates: LHA115- 
S 18 in the Small Magellanic Clo ud dClark et al.|2013 l 
and the Galactic object MWC 137 (Muratore et al.|2U15 ). 
For a better understanding of the evolution of these stars 
up to the supernova stage, it is hence essential to study 
B[e]SGs in great detail, and in particular their mass-loss 
behavior and the chemistry involved in the formation of 
the observed dense molecular and dusty rings or disk-like 
structures. 

The atmospheres of massive stars have an oxygen-rich 
composition, which means that the abundance of oxy¬ 
gen atoms greatly exceeds that of carbon atoms. The 
same trend is then expected in their disks, because they 
form from the material released from the stellar surface. 
Consequently, carbon atoms in the disks are locked in 
CO, which, as the most stable molecule with a dissocia¬ 
tion energy of 11.2 eV, forms as soon as the gas tempera¬ 
ture drops below the dissociation temperature of CO, for 
which a value of ^ 5000 K is typically quoted. The ex¬ 
cess oxygen atoms will form other molecules, and a rich 
chemistry is expected at larger distances from the star. 

The SiO molecule is the second most stable oxy¬ 
gen compound. Its binding energy of 8.0 eV is higher 
than that of other oxides or oxygen compounds such 
as NO, OH, or water, but lower than that of CO, and 
thus SiO should form in slightly cooler environments 
(TsiO < Tqo < 5000K). As with CO, the rather high gas 
temperature of the SiO forming region guarantees that 
high vibrational levels will be excited, causing coupled 
rotational-vibrational transitions, detectable via promi¬ 
nent band and band head features. And in fact, SiO 
bands are commonly observed in absorption in the at- 
mosph eres of Mira stars, cool giants, and red super- 


2 :iants dHinkle et al. 1976 Geballe et al. 1979 Aringer 
|et al.|[lM7f 1999). lii these objects, SiO is an important 
molecule, based on which the physical properties in the 
transition region between the outer atmosphere and the 
innermost cir cumstellar envelope can be studied (e.g., 
Ohnaka]|2014 ). 


'ib our knowledge, rotational-vibrational SiO bands in 
emission have only been reported for Supernova 1987A 
dAitk en et al. 1988 Roche et al. 1991 Liu & Dalgarno 
l994p. These emission bands appeared 161) days after 


the outburst and were not detectable beyond day 578. 


The rise of the SiO features had a time-delay of 48 days 
with respect to that of CO molecular bands. This is 
in agreement with the lower binding energy, requiring a 
cooler environment for the formation of SiO molecules by 
gas phase chemistry. The disappearance of the emission 
bands agrees with the onset of dust formation in the 
ejecta, which is interpreted as depletion of molecular SiO 
due to condensation of silicate dust, for which the SiO 
molecule is proposed to be a major building block within 
an oxygen-rich environment. 

Unlike the transient phenomena of a rapidly expand¬ 
ing, diluting ejecta of a supernova, the rings or disks 
around B[e]SGs seem to be quite stable structures of ac¬ 
cumulated material. This guarantees these disks as ideal 
chemical laboratories to study molecule formation and 
dust condensation. However, apart from CO band emis¬ 
sion, detected in most of these objects, and warm cir¬ 
cumstellar dust, no information exists about the chem¬ 
ical composition of these disks. Hence, it is necessary 
to begin searching for other molecular emission features 
that can be used to study the structure and kinematics 
of B[e]SG stars’ disks comprehensively. 

To pioneer this work, we selected a sample of Galactic 
B[e]SGs, with known CO band emission, to search for 
molecular rotational-vibrational band emission from SiO. 


2. OBSERVATIONS 

During the period of March-May 2014, we observed 
four Galactic B[e]SGs in the L-band wavelength region 
using the CRyogenic hi gh-resolution Infr aRed Echelle 
Spectrograph (CRIRES, Kaufl et al.||2004p on the ESO 
VLT UTl-Antu 8-m telescope. Details on the objects 
and their observations are given in Table[^ Data were 
obtained at four different standard settings with refer¬ 
ence wavelengths of 4048.6, 4060.1, 4135.5, and 4147.0 
nm, so that the total spectra cover a continuous wave¬ 
length range from 3.988-4.176/im. This range contains 
the first four band heads of the SiO first-overtone bands. 
We used the 0.4"slit, which provides a spectral resolution 
of i? ~ 50 000, corresponding to a velocity resolution of 
^ 6kms“^. 

Data reduction was performed with the ESO CRIRES 
pipeline (version 2.3.2). The observations were taken 
in an ABBA nod pattern for proper sky subtractions. 
Raw frames were corrected for bad pixels, flat fields, non¬ 
linearity, and then wavelength calibrated using OH lines 
in the standard star spectrum. A telluric standard star 
observation was similarly reduced. For telluric correc¬ 
tion, a B-type standard star was observed immediately 
after the target at a similar airmass. The IRAF^ task tel¬ 
luric was used to remove the atmospheric features. The 
final spectra were normalized and continuum subtracted. 


3. RESULTS 

We clearly detect emission from the first three band 
heads of the SiO first-overtone bands in CPD-52 9243 
(Fig.f^, and from the first band head in the spectra 
of the other three stars, CPD-57 2874, HD 327083, and 
HD 62623 (Fig.[|). 

^ IRAF is distributed by the National Optical Astronomy Ob¬ 
servatory, which is operated by the Association of Universities for 
Research in Astronomy (AURA) under cooperative agreement with 
the National Science Foundation. 
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In each star’s spectrum, the band heads display a blue- 
shifted shoulder and a red-shifted maximum with respect 
to the laboratory wavelength. Their separation is marked 
by the bar in Figs.[T] and Such structures are also 
typically observed in the band heads of CO, and are an 
unequivocal manifestation of rotational broadening. 

To conhrm that the observed features are emission 
from SiO molecules, we compute synthetic band head 
models. For this, w e modify our existing CO disk code 
(Kraus et al. 2000) by implementing the SiO li ne lists 
and F l mstein transition probabilities provided by |Barton| 
et al. (|2013 ). For each of our targets, information about 


the rotation velocity of the CO region and disk inclina¬ 
tion angles are listed in Table[^ The rotation velocities 
of the CO gas have been obtained from the broadening 
of the first band head of the CO first-overtone bands 
resolved in high-resolution spectra, and the disk inclina¬ 
tion angles have been determined from interferometric 
observations. From the wavelength separation between 
the position of the blue shoulder and the red peak of 
the first band head, we obtain an estimate of the SiO 
rotation velocity projected to the line of sight. This pa¬ 
rameter is refined during the fitting procedure. The fi¬ 
nal rotation velocities, corrected for disk inclination, are 
included in Tablej^ The high spectral resolution, com¬ 
bined with the sensitivity of the synthetic spectra to rota¬ 
tion velocity, guarantee that these values have a precision 
of ±lkms“^. Comparing the velocities obtained in the 
two molecular band emitting regions, we find a consis¬ 
tently lower value for SiO, suggesting that this molecule 
is formed within the Keplerian disk at larger distances 
from the star. 

For CPD-52 9243, where three observed band heads 
are available, we can also constrain temperature and col¬ 
umn density, for which we find Tsio = 2000 ± 200 K and 
-^SiO = (5 ± 2) X 10^^cm“^. With this column den¬ 
sity, most of the individual SiO rotation-vibration lines 
within the observed wavelength range are optically thick. 
Although the SiO rotation velocity in this object is only 
marginally lower, both temperature and column density 
are significantly l ower than the value s obtained from the 
CO modeling (see Cidale et al.|2012 1. This suggests that 
the two band head formation regions are physically dis¬ 
joint. 

For the other three objects, where solely the first band 
head is detected, it is only possible to constrain the rota¬ 
tion velocity, because different combinations of tempera¬ 
ture and density result in almost identical synthetic spec¬ 
tra. Hence, the hts shown in Fig.[^ are only to demon¬ 
strate that the observed features can indeed be assigned 
to SiO band emission, and to determine the SiO rotation 
velocity. 

Under the assumption of Keplerian rotation of the 
disks and using the stellar mass ranges of the objects 
(Table[^, the radii of both (CO and SiO) molecular 
rings can be computed. These values are included in 
Table|g Furthermore, considering that the circumstellar 
dust is mainly composed of silicates, we can calculate the 
dust evaporation distances. Assuming that the grains are 
spheres with a maximum size of ~ 1 /rm, and that the sil¬ 
icate evaporation temperature is ~ 1500 K, the minimum 
dust evaporation radii follow from the stellar luminosities 
(listed in Table[^ and are added to Table[^ Obviously, 
the much hotter molecular rings reside within the dusty 


disks. 


4. CONCLUSIONS 

We report on the detection of SiO first-overtone band 
emission in the high-resolution L-band spectra of four 
Galactic B[e]SGs, which are known to display strong CO 
band emission. While SiO bands are commonly seen in 
absorption from the atmospheres of cool stars and giants, 
they were detected in emission, so far, only in spectra 
of SN 1987 A during the post-explosion period from day 
160 to day 578. Hence, to our knowledge, this is the first 
discovery of SiO first-overtone band emission from hot, 
high-density environments of evolved massive stars. 

The high spectral resolution of our data allows us to 
precisely determine the rotation velocity of the SiO gas. 
For all of our target stars, this velocity is lower than that 
determined from the CO band forming region. Moreover, 
model fits to the observed three band heads in CPD-52 
9243 demonstrate that, in contrast to the CO band emit¬ 
ting region, the emission of the SiO bands originates from 
a cooler and less dense disk region, indicating that SiO 
molecules form at larger distances within the Keplerian 
disks. For all four objects, we find, in addition, that the 
SiO ring radii are smaller than the silicate dust evapora¬ 
tion distances. 

This finding is of great importance, because it demon¬ 
strates that molecules provide an excellent tool to study 
B[e]SG stars’ disks. The emission features of both 
molecules detected so far, CO and SiO, carry all the es¬ 
sential information about the physical properties (tem¬ 
perature, density) and kinematics of their formation re¬ 
gion. As different molecules require diverse conditions, a 
huge variety can be expected to form throughout the 
disk, filling the space between the hottest inner rim 
traced by CO bands, and the dust condensation region. 
Hence, it is essential and timely to search for emission 
features from many more molecules that are expected to 
form within the oxygen-rich environment of these disks. 
Moreover, molecules such as CO and SiO are also indis¬ 
pensable to study the disk structure at much larger dis¬ 
tances. With their numerous pure rotational transitions, 
these molecules are also suitable to trace much cooler en¬ 
vironments. The rotational transitions spread over the 
submm and radio regime, making the disks of B[e]SGs 
ideal targets to be observed with high spatial and spec¬ 
tral resolution facilities such as ALMA. Combining the 
results obtained from different molecular species tracing 
diverse temperature regions, and hence distances from 
the star, will greatly improve our knowledge and com¬ 
prehension of the disk structure, and help to understand 
the disks’ possible formation history. 

In addition, the discovery of SiO band emission has an 
even wider impact. Evolved massive stars, such as red su¬ 
pergiants, B[e]SGs, and supernovae (e.g., SN 1987A), are 
the main sources of silicate dust, for which SiO molecules 
are regarded as a major building block. Studying the 
physical properties of SiO gas in the disks of B[e]SGs 
might, hence, provide indispensable insight into the early 
stage of silicate dust formation within the oxygen-rich 
environment of these enigmatic stars. 
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Figure 1. Model fit (red) to the observed (black) three band heads of the SiO first-overtone band emission detected in CPD-52 9243. The 
arrows mark the wavelength of the unbroadened band heads, and the bar in the panel of the first band head marks the separation between 
the blue-shifted shoulder and the red-shifted maximum. 


Table 1 

Observed objects 


Object 

a (J2000) 

5 (J2000) 

K 

(mag) 

Obs. Date 

(K) 

(Mq) 

log L/Lq 

D 

(kpc) 

Ref. 

HD 62623 

07 43 48.469 

-28 57 17.37 

2.340 

2014-03-27 

8 500-9 500 

31-39 

4.8-5.2 

1.7±0.5 

1, 2 

CPD-57 2874 

10 15 21.971 

-57 51 42.71 

4.281 

2014-04-16 

9 000-10 000 

15-20 

6.00 

1.7±0.7 

3 

CPD-52 9243 

16 07 01.968 

-53 03 45.75 

4.443 

2014-03-22 

15 800 

17.4-18.6 

4.97 

3.44±0.8 

4, 5 

HD 327083 

17 15 15.374 

-40 20 06.74 

3.296 

2014-05-31 

11500 

60 

6.00 

4.80 

6 






20 000 

25 

5.00 

1.50±0.5 
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Table 2 

Rotation velocities and radii of the CO and SiO disk regions, and silicate dust evaporation 

distances 


Object 

i 

J^rot(CO) 

References 

■UrotCSiO) 

r(CO) 

r(SiO) 

revap (silicates) 


n 

(kms“^) 


(kms“^) 

(AU) 

(AU) 

(AU) 

CPD-52 9243 

46 

36 

1 

35.5 

11.9-12.7 

12.2-13.1 

16.1 

CPD-57 2874 

60 

130 

2, 3 

no 

0.8-1.0 

1.1-1.5 

52.7 

HD 327083 

50 

86 

4 

78 

3.0 / 7.2 

3.6 / 8.7 

16.7 / 52.7 

HD 62623 

38 

53 

3, 5 

48 

9.8-12.3 

11.9-15.0 

17.5 

References. 

— (1) |Cidale et al, 

!1(2012[; (2) 

IDomiciano de 

Souza et al 

ll2011|l; (3)[ 

Muratore et al.| 

; (4j Andruchow et al. (in preparation); (5) [Millour et al 
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Figure 2. Model fits (red) to the observed (black) first SiO band head. The arrow marks the wavelength of the unbroadened band head, 
and the bar marks the separation between the blue-shifted shoulder and the red-shifted maximum. 












































































































